Phase diagram of a spin ladder with cyclic four-spin exchange 
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We present the phase diagram of the 5 = 1/2 Heisenberg model on the two leg ladder with cyclic 
four spin exchange, determined by a combination of Exact Diagonalization and Density Matrix 
Renormalization Group techniques. We find six different phases and regimes: the rung singlet 
phase, a ferromagnetic phase, two symmetry broken phases with staggered dimers and staggered 
scalar chiralities, and a gapped region with dominant vector chirality or collinear spin correlations. 
We localize the phase transitions and investigate their nature. 

PACS numbers: 75.10.Jm, 75.40.Mg, 75.40.Cx 
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Frustrated interactions in quantum spin systems give 
rise to new and exotic phases but are little understood 
because of the inherent difficulties with competing inter- 
actions. A special type of frustration due to cyclic ex- 
change interactions was recently found to be important 
in the spin ladder material La ;E Cai4_ a ;Cu2404i [lj and in 
cuprate antiferromagnets such as La2Cu04 [2(. Similar 
multiple spin exchange interactions are known to be rel- 
evant for the nuclear magnetism of 3 He (3 and for the 
spin structure of a Wigner crystal 0] . 

In contrast to these systems on the triangular lat- 
tice, where already the nearest neighbor antiferromag- 
netic spin exchange is frustrating, in the bipartite two- 
dimensional (2D) square lattice model, relevant for the 
cuprates, the frustration only enters through the cyclic 
exchange term. It has been proposed that on the square 
lattice exotic magnetic phases with fractionalized excita- 
tions could exist 0. However, experience gained from 
quantum dimer models, showing a dimer liquid phase on 
a triangular lattice^ but not on the square lattice 
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FIG. 1: Phase diagram of the two leg ladder with cyclic 
four spin interaction. Squares denote first order phase tran- 
sitions, and the full circles indicate second order phase tran- 
sitions. The nature of the phase transition marked with the 
empty circle is presently unknown. The dashed line indicates 
a crossover line without a phase transition. 



indicates that on a bipartite lattice ordered phases are 
preferred and exotic spin liquid phases hard to realize. 

In our numerical studies of a cyclic exchange model we 
consider a ladder geometry, which is the simplest system 
on which cyclic exchange is possible. Consisting of two 
coupled chains, ladders already exhibit behavior remi- 
niscent of two-dimensional systems and can be investi- 
gated accurately using the density matrix renormaliza- 
tion group (DMRG) method 0]. While there exist a 
number of interesting ladder materials the aim of our 
study is not the discussion of these materials (which all 
seem to be in the Haldane phase). Rather, the results 
obtained give insight also into the properties of the 2D 
system. 

Previous numerical studies on ladders 0, 0, HH 
Il2^ were restricted to small cyclic exchange terms, only 
up to the critical point where the spin gap closes and a 
new phase appears. We go beyond this weakly frustrated 
regime and present the phase diagram at all strengths 
and signs of the bilinear and cyclic four spin exchange. 
Our SU(2) invariant Hamiltonian on a S = 1/2 two leg 
ladder is defined as follows: 

H = Ji J^Sfol) • 8(0!, 2) 

X 

+J|, ^S(i,»)-S(i + l,y) (1) 
+K J2 [P a {x,x + 1) + P^(x,x + 1)] , 

X 

where Jj_ (J\\) are the bilinear exchange constants on the 
rungs (along the legs) of the ladder and K denotes the 
coupling of the cyclic four spin permutation operators. 
This operator can be decomposed in terms of s pin oper- 
ators involving bilinear and biquadratic terms |l3| . We 
set Jj_ = J|| = J in the following and parameterize the 
couplings as J = cos(f?) and K = sin(f9). The energy 
scale \J J 2 + K 2 is set to one. 

The numerical algorithms we employ are exact diago- 
nalization (ED) of finite systems up to 2 x 16 sites with 
periodic boundary conditions (PBC), and the finite size 
version of the DMRG algorithm Q on systems with up 
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FIG. 2: (a) One of the two degenerate groundstates in the 
long range ordered staggered dimer phase. The spins are 
paired in singlets, (b) Spatial structure of the scalar chi- 
rality order parameter in one of the two symmetry breaking 
groundstates. The oriented triangle (the sign) denotes the 
orientation of the triple product (the sign of the order param- 
eter), (c) Correlation pattern in the dominant vector chiral- 
ity region. Oriented bonds stand for the cross product of two 
spins. Correlations between two such bonds are positive. 




FIG. 3: ED calculations of the staggered dimer structure fac- 
tor |17| (filled symbols) and the scalar chirality structure fac- 
tor (open symbols) in the ground state. Lines are a guide to 
the eye. The increase of the structure factor with system size 
is an indication for long range order. 



to 2 x 200 sites, keeping up to 1000 states and using 
appropriate open boundary conditions (OBC) 01 . We 
carefully checked the convergence of our results with re- 
spect to the lattice size and the number of states kept. 

The phase diagram in Fig. ^ summarizes our results. 
We now proceed to characterize the phases based on 
their order parameter or their slowest decaying corre- 
lation function (termed dominant correlation function) 
and then discuss the phase transitions, crossovers and 
universality classes. 

Rung Singlet Phase - We start the discussion of the 
phase diagram at 6 = 0, i.e the spin ladder with only 
antiferromagnetic bilinear couplings. Here the ground- 
state is unique and well approximated by the product of 
local rung singlets, hence the name rung singlet phase. 
All excitations are gapped and the correlation functions 
decay exponentially. The dominant correlations are the 
spin-spin correlations. A small positive K is sufficient to 
close the gap and to drive the system into a new phase 
H H IH [SI 113 We find that a negative K has a less 
pronounced effect. All correlation functions decay even 
faster and at 6 = -0.40 tt (J/\K\ = 0.30) we locate a 
first order transition to the ferromagnetic phase. 

Staggered Dimer Phase - Between 9 = 0.07 ± 0.01 7r 
(K/J = 0.23 ± 0.03) and 9 « 0.15 tt (K/J w 0.5) the 
system is in a valence bond crystal phase with a staggered 
dimer pattern and a finite gap to triplet excitations. The 
order parameter is: 

(S(:r -l,y). S(x, y) - S(x, y) ■ S(x + l,y)) (2) 

with a spatial structure as shown in Fig. [2] (a). This 
phase breaks the translation symmetry and has a twofold 
degenerate groundstate. Since this is a broken discrete 
symmetry, long range order (LRO) can exist even in one 
dimension. This order is seen in ED calculations (full 
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FIG. 4: Upper panel: local S(x,y) ■ S(x + l,y) expectation 
value on one of the two legs in the dimer LRO phase at 
8 — 0.12 n. The OBC render a direct measurement of this or- 
der parameter possible. Lower panel: long distance behavior 
of the scalar chirality correlations between equally oriented 
triangles in the scalar chiral LRO phase at 9 = 0.19 n. 



symbols in Fig. an d in DMRG resul ts ( upper panel of 
Fig. 0| . Since the boundary conditions |l5j in the DMRG 
calculations select one of the two possible groundstates, 
the order parameter can be measured directly and shows 
convincing evidence for long range staggered dimer order. 
Finally, the doubling of the unit cell in the symmetry 
broken phase manifests itself in degenerate singlets at 
momenta (0,0) and (n,ir) in the infinite system. We 
have confirmed the existence of two nearly degenerate 
states at these momenta in ED, with a small finite size 
splitting. 

Scalar Chirality Phase - For 9 larger than s» 0.15 7r 
the dimerization vanishes and we find a gapped phase 
with LRO in the staggered scalar chirality. The order 
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parameter is: 



(S(x,l).[S(x,2)xS(x + l,l)]) 



(3) 



and has a spatial modulation with wave vector (tt, n) (see 
Fig. [3 (b) for a pictorial representation) . This order pa- 
rameter breaks spatial symmetries and time reversal sym- 
metry, but not SU(2). LRO in this unexpected phase is 
seen as before in i) ED calculations of the corresponding 
structure factor (Fig. [31 open symbols), ii) DMRG cal- 
culations of the order parameter correlations, converging 
to a finite value at large distances (lower panel of Fig0J, 
hi) the existence of a (tt, 7r) singlet which is energetically 
close to the groundstate. The discrete symmetry break- 
ing in this phase suggests a finite triplet gap. We find 
that in both the dimerized and the scalar chirality LRO 
phases the triplet gap is finite but small [A(S — 1) < 0.1]. 

Dominant Vector Chirality Region - At 9 = 0.39 ± 
0.01 tt (K/J = 2.8 ± 0.3) we locate a second order phase 
transition to a short range ordered phase with a unique 
groundstate and a fully gapped excitation spectrum sim- 
ilar to the rung singlet phase. In contrast to the rung 
singlet phase the dominant groundstate correlations are 
not the spin-spin correlations, but correlations of the fol- 
lowing vector chirality order parameters: 

S(x,y) x S(x + l,y), S(x, y) x S(x,y + 1) (4) 

in staggered circulation arrangement [Fig. (c)]. This 
vector chirality is also called twist or helicity and can 
be regarded as a local spin current operator for bilinear 
Heisenberg Hamiltonians. Correlations are strong be- 
tween bonds on rungs or legs, but diagonal bonds are 
very weakly correlated. The vector chirality and the 
spin correlations shown in Fig. for DMRG calculations 
at 9 — 5tt/8 clearly demonstrate that the vector chi- 
rality correlations decay much slower (correlation length 
£ « 30) than the spin-spin correlations (£ « 12). The ex- 
istence of a small but finite gap in this region is confirmed 
by a finite size scaling of the triplet gap for the case of 
pure K (9 = vr/2) in the inset of Fig. (DMRG results). 
An extrapolation in l/L (L < 200) yields a lower bound 
on the infinite system gap: A(S ~ 1) > 0.016A". Compu- 
tations of dynamical spin and vector chirality structure 
factors in this region reveal another striking difference 
compared to the Rung Singlet phase: the lowest triplet 
excitation, with wavevector (tt, tt), is not a magnon (i.e. 
a spin flip excitation) but rather of the vector chirality 
type. It exhausts a large fraction of the spectral weight 
in the vector chirality structure factor. The presence of 
vector chirality correlations for pure K can be antici- 
pated at the classical level. The groundstate has a four 
sublattice structure where nearest neighbor spins are or- 
thogonal to each other ^^] and therefore maximize the 
vector chirality. Our results now suggest that the tran- 
sition from the classical to the 5 = 1/2 quantum case 
leads to short range order in both the spin and the vector- 
chirality correlations, with the latter becoming dominant. 
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FIG. 5: Semi- log plot of the long distance behavior of the 
staggered vector chirality correlations on the rungs and the 
spin-spin correlations for 6 = 57r/8. The spin correlations 
decay much faster than the vector chirality correlations. The 
open (filled) symbols for the spin correlations denote negative 
(positive) correlations. Inset: finite size scaling of the spin gap 
at 9 = 7r/2. 



At 9 0.85 7T (K/\J\ ss 0.5) there is a crossover region 
to dominant collinear spin correlations. 

Dominant Collinear Spin Region - In the proximity 
of the ferromagnetic phase boundary a short range or- 
dered region characterized by collinear (0, 7r) spin-spin 
correlations is observed, where spins on the same leg (on 
different legs) exhibit ferromagnetic (antifcrromagnetic) 
correlations. The system has a unique groundstate and 
a fully gapped spectrum. 

Ferromagnetic Phase - The last phase is the fully 
polarized ferromagnetic phase, located between two first 
order transitions at 9 — 0.94 tt [K/\J\ = 0.19) and 
8 = -0.40 tt (J/\K\ = 0.30). This phase extends be- 
yond the rigorous bounds — tt < 9 < —tt/2, inside which 
the ferromagnetic state minimizes the energy on each pla- 
quette separately. 

Phase Transitions and Universality Classes - We ap- 
ply the method of the Lieb-Schulz-Mattis twist operators 
to our system in order to precisely locate phases transi- 
tions and to discuss universality classes by considering 
the following quantity fl9| : 



(GS|exp 



x=l 



|GS), 



(5) 



where S z (x) = S z (x, 1) + S 2 (a;,2) and L is the system 
length. The infinite system values of this quantity con- 
ver ge t o ±1 or 0. In the Rung Singlet phase z L tends to 
+1 |ljj- The transition at 9 = 0.07 tt is signaled by a sign 
change in z L and can therefore be determined quite accu- 
rately. Our scenario of a transition from the Haldane uni- 
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FIG. 6: ED results of z L (definition in the text) as a function 
of 9. The infinite system values of this operator are ±1 or 0. 
For clarity not all system sizes are shown. 



versality class of the Rung Singlet phase to the dimerized 
phase is reminiscent of the Babudjian-Takhtajan (BT) 
type transition in the S=l bilinear-biquadratic chain |20j 
and this supports results obtained in 0, 0, 0] . The 
phase transition between the dimerized phase and the 
scalar chirality LRO phase (8 « 0.15 tt) has no visible 
signature in the behavior of z L . The rapid change of the 
groundstate correlations at this transition make it diffi- 
cult to discern between a first or second order phase tran- 
sition. The second order phase transition at 8 = 0.39 ir is 
detected again by a sign change in z L . The universality 
class of this phase transition remains for further study. 
There is no evidence for a phase transition between the 
dominant vector chirality and the dominant collincar spin 
regions but we see a smooth crossover with stable gaps 
instead. The ferromagnetic region finally shows vanish- 
ing z L and the transitions to it are of first order. 

Conclusions - The two leg ladder with cyclic four 
spin exchange reveals a very rich phase diagram. Be- 
sides determining the domain of stability of conventional 
phases such as the rung singlet phase, the collinear spin 
region and the ferromagnetic phase we have established a 
dimerized phase and chiral phases. Dimerized phases are 
common in frustrated spin chains but do not appear as 
generic phases of the diagonally frustrated two leg ladder 
[16} . However several authors have shown the existence of 
dimerized phases in ladder models including biquadratic 
terms [2U |22| and earlier studies of the present Hamilto- 
nian |8L licl Il2| conjectured a dimerized phase for large 
K/ J - which we have now confirmed. 

Unexpectedly we also found two chiral regions, one 
with long range order in the staggered scalar chirality 
and a second with dominant vector chirality correlations. 
Uniform scalar chirality phases have a long history in 
anyon superconductivity [23j and were discussed in the 
context of frustrated spin models [24| but were not seen 
or conjectured in the context of ladder models. 

Our results on ladders give insights also into the phase 
diag rams of cyclic exchange models on the square lattice 
|18| . since long range ordered phases on the ladder can 
be further stabilized when they are coupled to form 2D 



planes. We conjecture that the dominant vector chirality 
region leads to a T — long range ordered vector chiral- 
ity phase. This state can be regarded as a spin nematic 
which has long range order in the twist correlations, but 
no magnetic moment |2^| . The existence of such a phase 
has been conjectured, but a realization has not been ob- 
served in a microscopic model up to now. In fermionic 
terms this state is related to triplet (i-density wave states 
with staggered spin currents [26j. 

We would like to thank I. Affleck, M.P.A. Fisher, 
A. Kampf, C. Lhuillier, M. Sigrist, R.R.P. Singh, S. Todo, 
and S. Wessel for very fruitful discussions. G. S. and 
M. T. acknowledge support from Swiss National Science 
Foundation. We would like to thank S. R. White for help 
with the DMRG code. 



[1] 
[2] 
[3] 
[4] 



[5] 

[6] 

[7] 
[8] 
[9] 
[10] 

[11] 
[12] 
[13] 

[14] 
[15] 



[16] 
[17] 



[18] 
[19] 



[20] 
[21] 
[22] 



[23] 
[24] 
[25] 

[26] 



Matsuda et ai, Phys. Rev. B 62, 8903 (2000). 
Coldea et ai, Phys. Rev. Lett. 86, 5377 (2001). 
Misguich et al., Phys. Rev. Lett. 81, 1098 (1998). 
Chakravarty et ai, Phil. Mag. B 79, 859 (1999); 
Bernu et ai, Phys. Rev. Lett. 86, 870 (2001); 
Katano and D.S. Hirashima, Phys. Rev. B 62, 2573 

(2000) . 

T. Senthil and M. P. A. Fisher, Phys. Rev. B 63, 134521 

(2001) . 

R. Moessner and S.L. Sondhi, Phys. Rev. Lett. 86, 1881 

(2001) . 

S.A. Kivelson et al, Phys. Rev. B 35, 8865 (1987). 

S. Brehmer et al, Phys. Rev. B 60, 329 (1999). 

Y. Honda and T. Horiguchi, cond-mat/0106426 

K. Hijii and K. Nom ura, Phys. Rev. B 6 5, 104413 (2002). 

T. S. Nunner et ai, cond-mat/0203472 

M. Mueller et ai, Phys. Rev. B 66, 134423 (2002). 

The representation o f Pp +Pn 1 in terms of spin operators 

can be found in ref. |lq| for example. 

S. R. White, Phys. Rev. Lett. 69, 2863 (1992). 

We used standard open boundary conditions in all the 

regions but the dimerized and the scalar chirality LRO 

phase. In these cases we attached one more (5=1/2) spin 

at each end of the ladder in the same way as Ref. 16| . 

T. Hakobyan et al, Phys. Rev. B 63, 144433 (2001). 

We define the dimer structure factor in this context as the 

staggered sum of all leg dimer dimer correlations where 

the two bonds do not overlap. 

A. Chubukov et al, Phys. Rev. B 45, 7889 (1992). 

M. Nakamura and S. Todo, Phys. Rev. Lett. 89, 077204 

(2002) ; M. Nakamura and S. Todo, Prog. Theor. Phys. 
Suppl. 145, 217 (2002). 

L. A. Takhtajan, Phys. Lett. 87A, 479 (1982); H. M. 
Babudjian, Phys. Lett. 90A, 479 (1982). 
A. A. Nersesyan and A. M. Tsvelik, Phys. Rev. Lett. 78, 
3939 (1997). 

A. K. Kolezhuk and H.-J. Mikeska, Phys. Rev. Lett. 80, 
2709 (1998); M. J. Martins and B. Nienhuis, Phys. Rev. 
Lett. 85, 4956 (2000). 

X. G. Wen et al, Phys. Rev. B 39, 11413 (1989). 

D. Poilblanc et al, Phys. Rev. B 43, 10970 (1991). 

P. Chandra and P. Coleman, Phys. Rev. Lett. 66, 100 

(1991). 

C. Nayak, Phys. Rev. B 62, 4880 (2000). 



